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Abstract 

We show that the 2Aa signal of the leptonic flavor violating (LFV) Higgs boson decay h —>■ gr, as observed 
by the CMS collaboration recently, can be explained by a certain class of two-Higgs doublet models that 
allow controllable flavor-changing neutral current with minimal number of free parameters. We postulate 
that (i) the alignment limit is maintained, which means the lightest neutral scalar (h) has identical couplings 
to that of the Standard Model Higgs boson and (ii) the signal comes from two other neutral scalars, the 
CP-even H and the CP-odd A, almost degenerate with h at 125 GeV. We also show that (i) it is entirely 
possible that these scalars are hidden, apart from this LFV signal; (ii) the signal strengths of bb, 
and 77 around 125 GeV put severe constraints on the parameter space of such models; (hi) the constraint 
is further enhanced by the non-observation of processes like y —>■ ey, and we predict that the branching 
ratio of /i —>■ 67 cannot be even an order below the present experimental limit, highlighting the role it plays 
in forcing H and A to be near-degenerate; (iv) an enhancement in the production cross-section at 

around 125 GeV is expected in the gluon fusion channel, and should be observed during the next run of the 
LHG; (v) the branching ratio in the er channel is enhanced and is expected to be at least about 2%. The 
constrained parameter space and minimum number of free parameters, along with such strong predictions, 
make this model easily testable and falsifiable. 

PACS no.: 12.60.-i, 12.60.Fr, 14.80.Ec 


1 Introduction 


Leptonic flavor violating (LFV) processes do not take place in the Standard Model (SM); even with massive 
neutrinos, they are expected to be unobservably small, because the amplitudes are controlled by the tiny neutrino 
masses. Thus, observation of any LFV process is a smoking gun signal for physics beyond the SM [1]. No such 
signal has been observed so far in processes like £i —>■ £ 27 , t (-i —t ^ 2 ^ 1 , Mi -A M 2 iil 2 , where i and 

M stand for a generic lepton and meson, and the indices are arranged in such a way that the processes are both 
LFV and kinematically allowed. 

Recently, the CMS collaboration [2] found a 2.4(7 signal in the Higgs boson decay channel h -A p,T. This excess 
has been observed in both leptonic and hadronic final state channels of r; the weighted average is 

Br(h ^ far) = 0.841°;^^% , (1) 

and the upper limit at 95% confidence limit (CL) is Br(/i —)■ /rr) < 1.51%. 

There have been several attempts in the literature to explain this signal by introducing LFV couplings of the 
Higgs boson. This can be achieved by an extension of the scalar sector (and if necessary, the gauge and fermion 
sectors too), with some discrete (like S 4 [3] or A 4 [4]) or continuous (like gauged — Lr [5] or an 17(1)' with 
two scalar doublets [6]) symmetries, or supersymmetric Froggatt-Nielsen mechanism [7]. All these extensions 
necessarily introduce a number of new arbitrary parameters in the model. It was found [8] that an extension of 
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the scalar sector is imperative to explain the anomaly. As a concrete example, phenomenology of the type-III 
two-Higgs doublet model (2HDM) was considered in detail in [5,6,8-10], including some predictions for the 
model. The general feature of all these models is to predict at least one new scalar with LFV couplings that 
mixes with the SM doublet $ and the resultant mass eigenstate, which is dominantly d) with a small admixture 
of the LFV scalar, showing LFV signals while being in conformity with the SM for flavor-conserving decay 
channels. 

In this paper, we would like to explore the consequences of a particular class of 2HDM [11], first proposed 
by Branco, Grimus, and Lavoura (BGL) [12], which has tree-level flavor-changing neutral current (FCNC) 
interactions, appropriately suppressed by the entries of the quark or neutrino mixing matrix elements. The 
number of new parameters introduced in this class of models is minimal. In a certain limit (called the alignment 
limit) motivated by the LHG data, a particular type of the BGL model can not only explain the —>■ /rr signal, 
but also turns out to be extremely predictive. While a slight deviation from the alignment limit is still possible, 
resulting in the small-admixture explanation mentioned in the previous paragraph, there is another much more 
interesting possibility that we would like to explore. 

We speculate that in the alignment limit, the second CP-even neutral scalar H is almost degenerate with the 
Higgs boson h; so that the LFV signal comes from H and not h, which can decay in flavor-blind channels. But 
that is not the end of the story; H and the CP-odd neutral scalar A both contribute to /x —>■ ey (and other LFV 
processes as well, but this is most tightly constrained). If we want to keep Br(/x —>■ ey) within the experimental 
limit, H and A must also be near-degenerate. Thus, there are three neutral scalars sitting around 125 GeV, 
among which one is identical with the SM Higgs boson and the other two can have both flavor-conserving and 
flavor-violating couplings. We will discuss all the constraints on the parameter space, and show how one can 
successfully hide these two new scalars from the current LHC data. Concept of such degenerate Higgs bosons 
and their search strategies have been discussed in the literature [13,14]. One may note that signatures of light 
scalars of this model in the LHC experiments may not yet be observable [15]. 

If this model is true, we have a few tangible predictions. On the theoretical side, the ratio of two vacuum 
expectation values (VEV), commonly parametrized by tan/3, must lie in a narrow range, something like 0.4 < 
tan /3 < 2.8. This range can further be narrowed down with a more precise measurement of the neutrino mixing 
matrix. Another crucial input is the /i —>■ rr signal strength in the gluon fusion (ggF) production channel, and, 
as we will show, may point to an even narrower range of tan/3 centered around tan,d = 1. On the experimental 
side, the next generation experiments looking for ^ ey should be able to see it, as the rate must be close to 
the upper bound unless there is a fine tuned cancellation. Apart from that, there will be a significant excess in 
h —>■ TT (which, by default, also includes H ^ tt and A —)■ rr) in the gluon fusion channel, as the other scalars 
do not have any gauge couplings and hence cannot be produced in the vector boson fusion (VBF) channel. 
Precision study of the Higgs boson, either at the upgraded run of the LHC or at some future e+e“ collider, will 
certainly be able to test this model. For a current status of precision studies in the context of 2HDMs, we refer 
the reader to Ref. [16]. 


2 Formalism 


The scalar potential part of the BGL model [12,17] is like the other canonical 2HDMs with 17(1) symmetry, for 
whose formalism we will not go into detail. A soft breaking term is introduced in the scalar potential to prevent 
the appearance of a massless pseudoscalar [17]. For the Yukawa part, we will follow the notations of Ref. [12] 
as much as possible. The Yukawa part of the Lagrangian is given by 

2 3 

^y = -Y. Y. + (Apa;3^„(xa2$*)4 + + h.c.) , 

3 = 1 a,0=1 

( 2 ) 
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where $i and $2 are the two Higgs doublets, r'^(r^) and A'J(A^) are the Yukawa coupling matrices in the quark 
(lepton) sector, and m', d', e', u', Q' and L' stand for right-chiral up-type, down-type, charged lepton, neutrino,^ 
and the left-chiral SU(2) doublet quark and lepton fields respectively. The primes indicate that these fields are 
in the gauge basis; the mass basis fields will be denoted without primes. The generation indices are denoted by 
a and /3. Both and are 3x3 matrices. 

Now we impose the following BGL symmetry on the Lagrangian [12]: 

S '■ Q'k^ u'k ^ $2 —>■ (fc = either 1, 2 or 3). (3) 

All the fields except those which appear above, remain unaffected under S. Note that Eq. (3) violates lepton 
flavor universality by construction because we have singled out any one of the up-type quark fields and labeled it 
as k. Which up-type quark is labeled as k will lead to different models within the BGL class. Such a symmetry 
leads to FGNG in the down-quark sector. One can put the nontrivial transformation to d). instead of m). in Eq. 
(3), which leads to FGNG in the up-type sector. Since the FGNG constraints are much tighter for down-type 
quarks, the former class of BGL models are more predictive and we will focus on them only. We shall call our 
model U-, c- or t-type in accordance with fc =1, 2 or 3 respectively in Eq. (3). 

One can think of such an abelian symmetry for the lepton sector too. Note that k for leptons need not be the 
same as k for quarks. In fact, one can even have tree-level FGNG for up-type quarks and charged leptons. This 
is something that the reader should keep in mind when we discuss the possible constraints on these models. 

The GP-even^ neutral components of $1 and $2 mix with each other to give the mass eigenstates h and H, and 
the mixing angle is usually denoted by a. An intermediate basis {H°,R} (not the mass basis in general) can 
be obtained from the basis by a rotation through (3 = tan“^(z; 2 /ui) with the property that the state 

possesses exact SM-like couplings with the fermions and the vector bosons. The connection between these 
two bases is given by 

= cos{p — a)H + sin(/3 — a)h, R = — sin(/3 — a)H -\- cos{j5 — a)h. (4) 


Clearly, if we require that the lighter CP-even mass eigenstate, h, should posses SM-like couplings then we are 
led to the alignment limit^ sin(/I — a) « 1. This is the limit favored by the current 2HDM fits [18-25]. From 
now on, we will talk only about the states h, H, and A. 


The symmetry S of Eq. (3) leads to a very specific texture of the Yukawa matrices [12,26] which leads to the 
following Yukawa Lagrangian in the alignment limit: 


Cy = --h [dD^d + uDuu] - -H \d{NdPR + NlPL)d + u{NuPr + NIPl)u 

V V I 




d{NdPR - NlPL)d - u{NuPr - nIPl)v\ + I ^H+u {VN^Pr - NIvPl) d + h.c. I , (5) 


where we have suppressed the generation indices, and the unprimed fermion fields are in the mass basis. Dd 
and Du are diagonal matrices with the corresponding Yukawa couplings as the diagonal entries, so that the 
couplings of with the quarks are identical with the SM ones. The coupling matrices and Nd are given 
by 


Nu — diag{—TTiii cot j3 , rric tan /3 , m* tan /3} , 
Nu — diag{m„ tan /3 , —rric cot j3 , m* tan /3} , 
Nu = diag{m„ tan j3 , rric tan /3 , —rrit cot /3} , 


(iVd)T = tan/3 m^Sij - (tan/3 -I- cot/3)14* , 

{Nd)ij = tan/3 rmSij — (tan/3 -I- cot /3)V*iVcjmj , 
{Nd)ij = tan/3 rriiSij - (tan,d -|- cot /3)Vf*Vtjmj , 


( 6 ) 


^We introduce right-handed neutrinos and assume, for simplicity, that the neutrinos are pure Dirac particles. 

^We will assume all terms in the scalar potential to be real, so that the mass eigenstates are CP-eigenstates too. 

^If we require the heavier scalar, H, to be identical to then we are led to cos(/3 — a) ps 1. Since we will be assuming h and 
H to be quasi-degenerate, this limit is also a possibility. Whatever we comment here in the context of the limit sin(/5 — a) = 1 is 
equally valid also for the limit cos(/3 — q) = 1. More accurate measurements for the masses will be necessary to pinpoint the mass 
hierarchy and thereby to distinguish between these two limits. 
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where the indices i and j run over the down-type quarks d, s and &, and the superscripts indicate which type of 
model we are considering, i.e. which quark flavor has the nontrivial transformation. 


There are tight constraints coming from neutral meson mixing, which can in principle go through tree-level 
scalar exchange. For u- and c-type models, constraints from such mixings force mn = tua. For t-type models, 
such an exact degeneracy is not needed, in particular for low values of tan/3 [15, 17], but as we will see, a near 
degeneracy will be motivated from the LFV muon decays. 


A completely analogous formalism goes for the lepton sector. For massless neutrinos, one gets the coupling 
matrices by the following replacements: {Nu,Du) —0, F = 1, and Nd{Dd) —>■ Nf.{De), so that there is no 
leptonic FCNC. For massive neutrinos, one just replaces V by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
matrix U, and the (u, d) labels by {v, e). There can be three types of leptonic mixing models, which we will call 
^ 1 , 1 ^ 2 , and r'a type models. For example, in the z^i-type model, the coupling matrices are 


Ajy = diag{—TOj/j cot /3 , tan/3 , tan/3} , = tan jd rriidij — (tan/3 + cot P)U*^Uijmj . (7) 

The couplings of all scalars with charged leptons are given by Eq. (6) . The tree-level decay width for H (A) —>■ /rr 
can be written as 


T{H{A) fir) 




(lap -b \b\'^)imjj^A) - m 


m^) — 


- (w^ - (m^ 

« ^TOi/(A)(|ap + |5p), (8) 

where the interaction Lagrangian is generically written as 

Afint = Jl{aPL + bPii)TH{A) -b h.c. (9) 

The expressions for the Yukawa couplings a and b for H and A can be taken directly from Eq. (5): 

H iA 

= --M [{Ne)^rPR + {N,)%Pl] ^ + —M [{Ne)^rPR " {Ne)%PL] T + h.C. , ( 10 ) 

where the expressions for (fVe)/iT and can be obtained for i/i-model from Eq. (7), as 

{Ne)iir = - (tan/3 -b cot^) U^^Uirrur , (fVe)rp = “ (tan/3 + cot/3) U^^Uirm^ , (11) 

with Ui^ = — sin 6 >i 2 cos 023 — cos 0 i 2 sin023 sin 0 i 3 e“*‘^ and Uit = sin 0 i 2 sin 023 — cos 0 i 2 cos 023 sin 0 i 3 e“®'^. The 
angles and phase are the same as used in the standard parametrization of the PMNS matrix. 


A very tight constraint on the parameter space comes from the radiative decay /x —>■ eq. The amplitude for the 
LFV decay £i —>■ £jj can be written as 

= -^rndj{ia^,q''){CLPL + . ( 12 ) 

Using this, we may write the expression for the BR as [27] 

BR(£,^^,7) = —(jCip + jCxxp) . (13) 

Since the charged scalar loop contribution depends on the ratio (m^./m^^.), we neglect it on account of tiny 
neutrino masses. The dominant contributions to Cl.r will come from the neutral scalar loops mediated by Pd 
and A. Detailed expressions for the contributions from neutral scalar loops, for the process b —> sy, already 
appear in the Appendix of Ref. [15]. From these, corresponding expressions for /i —>■ ey can be easily obtained 
by straightforward replacements of the CKM elements by the appropriate PMNS elements and the down-type 
masses by the corresponding masses of the charged leptons. 
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BRjfj -> ey) < 5.7x 10 
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Figure 1: (L) The allowed region in the tan P^niA plane coming from the non-observation o/^ 67 . The 

horizontal width comes from an almost fine-tuned cancellation between H and A. (R) The Br(fj, —>■ ej) shows 
more clearly the allowed region and the cancellation. The shaded region is still experimentally allowed. 


3 Analysis 


For our analysis, we will assume the BGL model to be t-type in the quark sector (u- and c-type models are very 
tightly constrained from flavor data [15]) and vi type in the lepton sector, and will call it ti^i model for brevity 
(we will use the same type of nomenclature for other models in the BGL class). Let us first try to justify our 
choice. 



Log,o(tan/?) 

Figure 2: The combined branching ratio for H,A ^ pLT in the vi type model. The lower shaded region shows 
the 95% upper limit for the branching ratio. In our evaluation of the BR, we have used F®'^ =4.1 MeV as the 
denominator to consistently compare our results with the bounds given in [2J. 

Our starting point is the assumption that the alignment limit is exact, as far as experimental precision goes. 
In this limit, h has only flavor-conserving couplings, and these couplings (also the gauge couplings) are exactly 
the same as those in the SM. Thus, the /rr final state, apparently originating from a 125 GeV scalar, cannot 
come from h. We postulate that the fir signal is the first hint of a quasi-degenerate nonstandard Higgs (in this 
case H with mn = 126 GeV). Under this assumption, we have checked that in the V 2 and 1^3 models, the decay 
H fXT overshoots the experimental data for any values of tan /?, while for the vi model, the decay is found to 
be under control for a narrow range of tan/?. As we will discuss later, the experimental bound on BR(/i —> ey) 
compels the pseudoscalar (A) to be equally light when mn = 126 GeV is assumed for the ni model. Thus, we 
are forced to think of a scenario where all three neutral scalars of the 2HDM are quasi-degenerate at about 125 
GeV. Therefore, it is important to include the contributions from both H and A while comparing with the LHC 
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( 14 ) 


signal strengths into different decay channels. To this end, we define all the /r-parameters as 

<j{pp h, H, A ^ F) 
a{pp hsM -t F) 

for any final state F accessible from the decay of the Higgs boson. Of course, F[ and A can contribute to flavor- 
conserving fermionic final states, or even to 77 . However, in the alignment limit, they do not have couplings 
of the form H{A)VV, where V = IT, Z is a vector boson. So there will neither be any production of i? or H 
through vector boson fusion, nor any decay to WW* or ZZ*. 

Now we discuss the constraints in greater details. With leptonic flavor mixing, both H and A can contribute to 
LFV processes. The tightest constraint comes from p, —>■ ey, which is shown in Fig. 1. This shows that, barring 
an unnatural cancellation between iJ and A contributions, the data forces A to be near degenerate with h and 
F[. We have also checked that, for the allowed region from p ey, the BRs of the processes r —?> ey and 
T ^ p'y are always < 10“^^ which is four orders of magnitudes below the current experimental limits [28]. So, 
we have the tantalizing possibility that there may actually be three states within the 126 GeV resonance. 

As far as the decay /r —> 3e is concerned, we have checked explicitly that over the entire parameter space that 
we talk about, the branching ratio at tree-level is at most of the order of 10 “^^, with almost degenerate scalars. 
This is about 5 orders of magnitude smaller than the current experimental upper limit of 1.0 x 10“^^ at 90% 
CL. The loop-level decay is suppressed compared to /r —>■ ey by at least one more power of a [29], so this process 
is well under control. 

Once we accept H and A to be degenerate with h, the pr signal has to come from both these scalars In Fig. 2 
we show the branching ratio for this channel in the i>i type model, with mn = 126 GeV and ruA = 124 GeV. 
The width of the parabola-like region comes from the 3cr experimental range of the PMNS angles [31] together 
with the variation of the Dirac CP phase, S in the range [0, 27r]; and a more precise determination of them will 
help in thinning out the parabola. The nature of the plot is not surprising because, as can be seen from Eq. 
( 6 ), the FCNC couplings come with a prefactor of (tan /3 -|- cot /3) which has a minimum at tan /3 = 1. What the 
figure shows is that tan/3 has to be extremely constrained, 0.37 < tan/3 < 2.8. For other i/i type models, the 
parabolic shape is maintained but the lowest point of the parabola lies well above the experimental range (this 
can be attributed to relatively larger magnitudes of the PMNS elements which are involved in the tree-level 
pT couplings for ^' 2,3 models), which forces us to consider only the i^i type model. It is worth noting that our 
scenario predicts BR{X —)■ pr) > 0.5% which can be tested in the next run of the LHC, where X denotes the 
summed contributions from H and A decays. 

Now that we have laid out the constraints, it is time to see why H and A could have been missed at the LHC. 
The only observable final states where one can have signals for these scalars are bb, yy, and rr, because these 
scalars do not have any trilinear gauge couplings. The experimental numbers, taken from Refs. [32-36], are, at 
the 1(7 level, 

= 1.93 ±1.45 (ATLAS), 0.52 ± 0.9 (CMS), 
pf^ = 1.32 ±0.38 (ATLAS), 1.12 ±0.37 (CMS), (15) 

while there are no significant bounds for p^^^ ■ 

The /r-values for these three channels, as a function of tan/3, are shown in Fig. 3; also shown are the 3cr bands 
for the ATLAS and CMS experiments. As the model is a t-type in the quark sector, bb or even yy constraints 
are relatively relaxed for tan/3 > 1. In fact, there is no noteworthy measurement of the signal strength in the 
bb channel with ggF tagging because of the huge background. The rr channel, on the other hand, shows a 
significant enhancement. Note that because of the nature of the PMNS coupling, the branching ratio depends 
on the leptonic CP-violating phase S', the least enhancement is expected for 5 = ir. 

Thus, apart from the rr channel, there is no reason why the new scalars should have been already seen at 
the LHC, if tan/3 > 1. Unfortunately, we cannot go to large values of tan/3 {e.g. 10) because that will be in 

^For exactly degenerate states, the decay widths may change because of the interference effect [30]. This, however, is a fine-tuned 
possibility which we will not enter into. 
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Figure 3: ^-values (for the definition, see Eq. (14)) for different final states in the ggF production channel, as a 
function of tan j3 for tvi-model. Contributions from all the three neutral scalars (h, H and A) have been included 
in the estimation of the signal strengths. For the plots we have assumed = 125 GeV, mn = 126 GeV and 
ruA = 124 GeV, and the central values for various PMNS matrix elements. Also shown are the ia limits on the 
fi-parameters from the ATLAS and CMS experiments. There are no significant numbers for h ^ bb in the ggF 
channel. 


conflict with H{A) pLT data. The ATLAS data on ditau signal strength in the ggF channel is, till now, pretty 
inconclusive [36], and supports as large as 5. The CMS data [33] is much more precise (although the ggF 
and the VBF channels are not always differentiated very successfully [37]). Anyway, all the relevant /x-values 
are well within the experimental 3 ct range; if we take the range at the 2a level, ~ 2.3, which forces the 
parameter space to a region near tan/3 = 1 and large nonzero 5. We find /x^gF ~ 5 fQj. tan/3 « 1, which is still 

within the 2a allowed range of both ATLAS and CMS [33,36]. This is another prediction which can be tested 
in precision Higgs studies. 

In passing, we note that while finding a similar excess in the e/x channel will be tough due to the smallness of 
the masses of the leptons involved, our scenario does predict a branching ratio for the er channel to be above 
2% (see Fig. 4), which can be used, among other things, to falsify the model. 
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Figure 4: Prediction for the branching ratio in the 
er channel. The PMNS angles have been varied 
within their 3a experimental range. 


Figure 5: Constraint on at 99% G.L. from the 
experimental limit on AT. 


The charged scalar merits a brief discussion. With nearly degenerate H and A, the strongest indirect bound 
on the charged scalar mass ( 171 ^+), for the BGL model at hand, comes from the oblique T-parameter. In the 
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2HDM alignment limit, the new physics contribution to the T-parameter can be expressed as [38,39] 


AT 

with, T{x,y) 


•——: 2 , TO^) + -F(m^+, m^) - 

iDTrsin OwM^ 

I 0 ior X = y . 


ttA)] . 


Taking the new physics contribution to the T-parameter as [40] 


AT = 0.05 ±0.12, 


( 16 ) 

(17) 


(18) 


we show the 3cr range for the allowed charged scalar mass in this model in Fig. 5. We find that the value of 
mH+ can rise as high as 250 GeV. Since the LHC is most sensitive for mH+ up to 160 GeV [41,42], a charged 
scalar heavier than that can remain well hidden from the current reaches of LHC. 


4 Summary 


We have shown that the recently observed excess in the yr channel coming out of the 126 GeV resonance can 
be interpreted in a radically different way than the canonical explanation of having one more scalar with LFV 
couplings which mixes with the SM doublet and the one of the resultant CP-even neutral mass eigenstates 
showing the LFV decay signal due to the slight admixture. We took a 2HDM with nonzero FCNC, and note 
that this model can also be tuned to the canonical explanation mentioned above to generate the LFV decays, 
with hardly any originality. Taking a different approach, we assumed that the alignment limit is exact, and 
explore the possible consequences. This is how our logic went. 


• In the alignment limit, the SM Higgs boson h cannot have any tree-level FCNC coupling. The yr signal, 
therefore, must be coming from another scalar, which is almost degenerate with h. The splitting has to 
be less than the resolution for two nearby mass peaks. 

• The 2HDM of Branco, Grimus, and Lavoura has at least two neutral scalars with tree-level FCNC, one 
being the CP-even H and the other, CP-odd A. We, therefore, assume H to be nearly degenerate with h. 

• This poses a serious problem with the non-observation of /r —>■ ey. This process receives contribution from 
both H and A mediated diagrams, but the amplitudes come with opposite signs. The way out is to take 
A to be nearly degenerate with H, so that the contributions more or less cancel. 

• H and A cannot be differentiated at the LHC. Thus, we have to take both H ^ jir and A ^ yr into 
consideration while calculating the branching ratio. Due to the nature of the coupling, the combined 
branching ratio has a minimum at tan /3 = 1 and grows on either side. For r '2 and type models, the 
lowest possible value of the branching ratio is way above the data, so we are forced to consider the vi 
model. This gives an allowed range of tan /? as 0.4 < tan /3 < 2.8. If we take a t-type model, all constraints 
coming from hadronic physics are satisfied in this range of tan/?. 

• H and A do not have any trilinear gauge couplings in the alignment limit. Therefore, the only way to 
produce them at the LHC is through gluon fusion (they can be radiated off a top, though). In the gluon 
fusion channel, there is no significant constraint for bb final states. The 77 final states are significantly 
enhanced by top loops for tan/? < 1 but is well under control for tan/? > 1. The only important channel 
is rr, which receives a significant enhancement. However, with the present state of the data, signatures 
for the new scalars in these channels can still be below any statistical significance. 

• Most importantly, the model is testable and falsifiable through various means. First, a more precise 
determination of the branching ratio in the yr channel can rule out this model, if the upper limit drops 



below the lowest possible value of the branching ratio. Note that the model has no free parameters once 
one specifies tan (3 (the neutral scalars are all nearly degenerate and the charged scalar is not relevant here; 
however, with mH ~ vriA in the alignment limit, cannot be too far away from mt for the theory to be 
consistent with the electroweak precision observables, in particular the T-parameter ® [17]) and therefore 
no way to play with the minimum branching ratio. Second, the model predicts a large branching ratio 
for /i — >■ ey, almost close to the experimental upper limit, which should be seen in the next generation 
of experiments, unless there is an unnatural cancellation between two competing amplitudes. Third, the 
model predicts a large branching ratio for the er final state, whose lowest value is about 2%. Fourth, and 
the most clinching, is a precise determination of h —>■ rr in the gluon fusion channel. If there is hardly 
any scope for a large enhancement, this model is ruled out. 
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